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ABSTRACT 



We present mid-infrared (8-33 pm) observations of the binary L5- Trojan sys- 
tem (617) Patroclus-Menoetius before, during, and after two shadowing events, 
using the Infrared Spectrograph (IRS) on board the Spitzer Space Telescope. 
For the first time, we effectively observe changes in asteroid surface tempera- 
ture in real time, allowing the thermal inertia to be determined very directly. 
A new detailed binary thermophysical model is presented which accounts for 
the system's known mutual orbit, arbitrary component shapes, and thermal 
conduction in the presence of eclipses. 

We obtain two local thermal-inertia values, representative of the respective 
shadowed areas: 21 ± 14 J s'^/^R-^m-^ and 6.4 ± 1.6 J s-^^K'^m-^. The aver- 
age thermal inertia is estimated to be 20 ± 15 J s"-'^/^K"^m"^, potentially with 
significant surface heterogeneity. This first thermal-inertia measurement for a 
Trojan asteroid indicates a surface covered in fine regolith. Independently, we 
establish the presence of fine-grained (< a few pm) silicates on the surface, 
based on emissivity features near 10 and 20 pm similar to those previously 
found on other Trojans. 

We also report y-band observations and report a lightcurve with complete 
rotational coverage. The lightcurve has a low amplitude of 0.070 ± 0.005 mag 
peak-to-peak, implying a roughly spherical shape for both components, and 
is single-periodic with a period (103.02 ± 0.40 h) equal to the period of the 
mutual orbit, indicating that the system is fully synchronized. 

The diameters of Patroclus and Menoetius are 106 ±11 and 98 ± 10 km, re- 
spectively, in agreement with previous findings. Taken together with the sys- 
tem's known total mass, this implies a bulk mass density of 1.08 ±0.33 g cm"^. 
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significantly below the mass density of L4-Trojan asteroid (624) Hektor and 
suggesting a bulk composition dominated by water ice. 

All known physical properties of Patroclus, arguably the best studied Trojan 
asteroid, are consistent with those expected in icy objects with devolatilized 
surface (extinct comets), consistent with what might be implied by recent 
dynamical modeling in the framework of the Nice Model. 



Keywords: Asteroids, Composition; Asteroids, Surfaces; Eclipses; Infrared 

Observations; Trojan Asteroids 
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1 Introduction 



1.1 Thermal inertia 



Thermal inertia is a measure of the resistance to changes in surface tempera- 
ture. A hypothetical zero-thermal-inertia asteroid would be in instantaneous 
thermal equilibrium with insolation and display a prominent sub-solar tem- 
perature peak. In general, thermal inertia causes the surface temperature dis- 
tribution to have reduced contrast and be asymmetric with respect to the 
sub-solar point, with the maximum shifted to the afternoon side. The influ- 
ence of thermal inertia on the temperature distribution is maximum when the 
spin vector is perpendicular to the solar direction and zero if it points towards 
the Sun. 



Thermal inertia is defined as F = y/K,pc, with thermal conductivity k, surface 
bulk mass density p, and heat capacity c. The range in which p and c can plau- 
sibly vary is much more restricted than for n, which can vary by several orders 



of magnitude (see, e.g., Mueller, 2007, Sect. 3.2.2.b for a detailed discussion). 



Examples of low and high thermal inertia cases would be a dusty, thermally 
insulating surface such as that of the Moon, and a rocky, thermally conduc- 
tive surface, respectively. It is qualitatively clear that for fine-grained regolith 
(where the typical grain size is not large compared to the penetration depth 
of the heat wave), thermal conductivity and thermal inertia decrease with de- 
creasing grain size. This is confirmed by laboratory studies under a simulated 



Martian atmosphere (Presley and Christensen, 1997), but no quantitative data 



are available for vacuum conditions. 
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The thermal inertia of the largest few main-belt asteroids is known to be 



low and indicative of a surface covered in fine regolith (e.g. Spencer, 1990 



Miiller and Lagerros, 1998). It was recently found that the thermal inertia of 



D < 100 km main-belt asteroids and near-Earth asteroids is increased relative 
to large main-belt asteroids, more so for smaller objects, but that it stays sig- 



nificantly below values expected for bare- rock surfaces (Mueller, 2007 Delbo 



et al. , 2007 Delbo' and Tanga, 2009). No quantitative measurements of the 



thermal inertia of Trojans have been published. 

The thermal inertia of asteroids is typically determined indirectly, from multi- 
filter observations in the challenging mid-infrared wavelength range at widely- 
spaced phase angles, such that the sub-observer point is located at different 
rotational phases ( "local times" ) and information on the diurnal temperature 
distribution can be obtained. Depending on the apparent target motion, this 
may require observations stretched out over a significant time span. 



Eclipse-induced shadowing, on the other hand, offers an opportunity to deter- 
mine the thermal inertia more directly, by monitoring the shadowing-induced 



cooling and the later warming up. This technique was pioneered by Pettit 



and Nicholson (1930), who determined the thermal inertia of the Moon for 



the first time from thermal observations during a total lunar eclipse. Morri- 



son and Cruikshank (1973) report observations of the Galilean satellites as 



they were eclipsed by Jupiter; Neugebauer et al. (2005) report observations 



of lapetus during eclipses by Saturn's rings; eclipse observations of Saturn's 



satellites using the Cassini spacecraft are ongoing ( Pearl et al. , 2008 ) . 



Asteroids are not frequently shadowed by planets. Binary asteroids, however, 
can undergo mutual eclipses, during which one component shadows the other. 



7 



Recent progress now allows mutual eclipse events in an increasing number of 
binary asteroid systems to be reliably predicted. 



1.2 Binary Trojan asteroid Patroclus 



Our target, (617) Patroclus, is one of the few known binary systems in the 
population of Trojan asteroids, which are in 1:1 orbital resonance with Jupiter, 
close to the Lagrangian points L4 and L5. Trojan orbits are stable over most 



of the age of the Solar System (Levison et al. , 1997). Their origin is currently 



under debate: While they were long believed to have formed near their present 



position (see, e.g., Marzari et al. 2002), Morbidelli et al. (2005 2009) argue in 



the context of the 'Nice model' that Trojans were captured by Jupiter during 
a chaotic phase in the early Solar System, which was caused by a mean- 
motion resonance between Jupiter and Saturn. In this model, Trojans share 
a volatile-rich parent population with comets and Trans-Neptunian Objects 
but have spent parts of their lifetime at relatively small heliocentric distances 
(potentially enough to leave their surface devolatilized) before being trapped 
in their current orbits. 



Trojans have generally low albedos of p\ ~ 0.04 and virtually featureless, 
highly reddened reflection spectra in the visible and near-IR wavelength ranges; 



in both respects, they resemble cometary nuclei (see Dotto et al. , 2008, for 



a recent review). Mid-IR spectra of three Trojans have revealed the presence 



of fine-grained silicates on the surfaces (Emery et al. , 2006). Trojans smaller 



than some 70 km in diameter appear to be collisional fragments, while larger 
bodies appear to be primordial "accretion survivors," i.e. bodies whose current 
form and internal structure have remained unchanged since the time of their 
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formation (Binzel and Sauter 1992 Jewitt et al. , 2000). 



(617) Patroclus was discovered in October 1906 by August Kopff. Following 
(588) Achilles, which had been discovered in February 1906, Patroclus was 
the second known Trojan and the first known object in the L5 swarm. With 
an absolute magnitude in the HG system (V magnitude normalized to a stan- 



dard geometry; see Bowell et al. , 1989) of if = 8.19, Patroclus is among the 



largest Trojans. Based on IRAS observations, Tedesco et al. (2002) report 



a diameter of D = 140.9 ± 4.7 km, corresponding to a geometric albedo of 



Py = 0.047 ±0.003 (see eqn. 1 below). Fernandez et al. (2003) obtained Keck 
12.5 pm observations of Patroclus and derived a diameter of 166.0 ± 4.8 km 
{py = 0.034 ± 0.004). Note that the quoted diameters equivalent 
diameters of the system as a whole. is related to the components' 
diameters, Di and D2, through D_4^ = Di - + D2^. 



Patroclus was found to be binary by Merline et al. (2001). The system's mu- 



tual orbit was first determined by Marchis et al. (2006a) based on spatially 



resolved adaptive-optics observations. Subsequently the secondary component 



was christened Menoetius. Berthier et al. (2007) provide an updated orbit 



model (consistent with that by Marchis et al. ) based on additional spatially 
non-resolved optical photometry measurements during mutual events. The 



Berthier et al.| model is used throughout this paper. The system's mutual or- 
bit is purely Keplerian (no precession) and circular (eccentricity < 0.001). 
The center-to-center separation of the two components is 654 ± 36 km, corre- 
sponding to a maximum angular separation around 0.2". The J2000 eclip- 
tic coordinates of the spin-pole orientation are A = 241.37 ± 0.33°, (3 = 
—57.35 ± 0.36°, the orbital period equals 102.94 ± 0.11 h. This orbit model 
reproduces the adaptive-optics observations and the morphology of eclipse- 
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induced lightcurves very well. Predicted and measured eclipse timing may 
differ by up to a few hours (= a few percent of the orbital period). Work is 
under way to resolve these slight timing differences. No discernible change in 
the model's spin-axis orientation is expected to result from this. 



Marchis et al. (2006a) report a component-to-component magnitude difference 



of ~ 0.17 mag; the same magnitude difference was found for two filters cen- 



tered at 1.6 and 2.2 pm, respectively, which [Marchis et al.| take as an indication 
of similar surface composition. Assuming identical albedo, the corresponding 
diameter ratio is D1/D2 ~ 1.082. 



Through Newton's parametrization of Kepler's Third Law, the orbital param- 



eters imply a total system mass of (1.20 ± 0.20) x 10^* kg. Marchis et al. used 



the size estimate by Fernandez et al. (2003) to infer a bulk mass density of 
only 0.8^0;^ g cm"^, compatible with a composition dominated by water ice and 



moderate porosity. Berthier et al. determined the components' diameters inde- 



pendently, from eclipse timing, to be Di = 112 ±16 km and D2 = 103 ±15 km, 
and a bulk mass density of 0.90l|]i^ g cm . 



1.3 Work presented in this paper 



The Marchis et al. (2006a) orbit model allowed a series of mutual events to 
be predicted, including two shadowing events in June 2006 which we observed 
using the Spitzer Space Telescope. We observed the system's thermal emission 
before, during, and after the events (see sect. 3). Although the system was not 
spatially resolved, Spitzer's high sensitivity allowed us to effectively observe 
the eclipse-induced cooling and later warming up of shadowed surface elements 



10 



8.2 



8.3 - 



s 



8.4 - 





^ 8.5 - 



8.6 



□ 1996 Apr 9.2 
o 1996 Apr 10.2 
A 1996 Apr 11.3 
+ 1996 Apr 12.3 



X 1996 Apr 13.2 
1996 Apr 15.2 
* 1996 Apr 16.2 
■ 1996 Apr 17.3 



• 1996 Apr 18.1 
A 1996 Apr 20.2 
P,^ = 103.02 hr 
T„ = 50186.0 MJD 



8.7 



0.0 



0.2 



0.4 



0.6 



0.8 



1.0 



1.2 



Rotational Phase 



Fig. 1. Composite visible lightcurve from our April 1996 observations. Data points 
at rotational phases 0.0-0.2 are repeated at 1.0-1.2 for clarity. Synthetic lightcurves 
are superimposed, from a 2nd order Fourier fit to the data (solid line) and from a 
1.07:1.0:1.0 ellipsoid model (dash-dotted line). See text for details. 

in real time. No such observations have been performed before. We present 
the first direct measurement of the thermal inertia of an asteroid. A new 
thermophysical model for eclipsing binary asteroids has been developed to 
facilitate the data analysis (see sect. 5). 

We also obtained Patroclus' mid-IR emissivity spectrum from our Spitzer 
data (see sect. 4), and its visible lightcurve from ground-based observations 
(see sect. 2). 

2 Visible lightcurve observations 

We observed (617) Patroclus using the 61 cm Bochum Telescope at the Euro- 
pean Southern Observatory in La Silla during 10 nights in April 1996; see table 

1 for observational details. The observations were performed in the Johnson 
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Table 1 

Optical observations: Observation mid time, observer-centric position (A and /? in 
J2000 ecliptic coordinates), solar phase angle a, heliocentric distance r, and geo- 
centric distance A. V{l,a) denotes the observed mean magnitude corrected for 
heliocentric and geocentric distance, but not for phase angle. 



UT date 


A 


P 


a 


r 


A 


V{l,a) 


(Apr 1996) 


(deg) 


(deg) 


(deg) 


(AU) 


(AU) 




9.2 


212.5 


6.8 


2.4 


5.7592 


4.7839 


8.49 


10.2 


212.4 


6.8 


2.3 


5.7585 


4.7795 


8.49 


11.3 


212.3 


6.7 


2.2 


5.7578 


4.7795 


8.49 


12.3 


212.1 


6.7 


2.0 


5.7572 


4.7715 


8.48 


13.2 


212.0 


6.7 


1.9 


5.7566 


4.7682 


8.45 


15.2 


211.8 


6.6 


1.6 


5.7553 


4.7621 


8.44 


16.2 


211.6 


6.6 


1.5 


5.7546 


4.7596 


8.44 


17.3 


211.5 


6.6 


1.3 


5.7539 


4.7573 


8.42 


18.1 


211.4 


6.5 


1.2 


5.7533 


4.7555 


8.42 


20.2 


211.1 


6.5 


1.1 


5.7520 


4.7525 


8.43 



V band with the DLR MKII camera, using a Tektronix Ik x Ik CCD. 

Differential photometry was performed against field stars, which were subse- 
quently absolutely calibrated against observations of standard stars present in 
near-by fields. Fluxes were extracted through synthetic aperture photometry. 
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using the AstPhot program (Mottola et al. , 1995). The typical la uncertainties 



in the differential photometry were well below 0.01 mag; the typical absolute 
photometric errors were of the order of 0.02 mag RMS. 

The synodic rotation period was determined by using the Fourier analysis 



technique described by Harris and Lupishko (1989). This technique seeks the 



period that, in the least-squares sense, best fits the observations. In order to 
remove the effects of changing brightness on subsequent nights due to solar 
phase effects, or to uncertainties in the nightly absolute photometric calibra- 
tions, the method enables grouping of observations (with one group usually 
comprising observations from one night) and estimates magnitude shifts for 
each group that minimize the residuals. In our case it became soon apparent 
that the time spanned by each night of observation was much shorter than 
the rotation period of Patroclus. As a consequence, when folded into a com- 
posite, single-night data stretches would only sporadically overlap with each 
other, thereby providing little constraint for the solution of the single-night 
magnitude shifts. For this reason, whenever possible, we grouped the obser- 
vations in sets of two consecutive nights, which resulted in a total of 6 groups 
with a good mutual overlap. Subsequently, the magnitudes within each data 
set were reduced to a reference phase angle (chosen as the phase angle of the 
observations closest to the average phase angle in each data subset) by using 



the HG model (Bowell et al. , 1989) and a G value of 0.12. Because the phase 
angle change rate at the time of the observations was only of the order of 0.1° 
per day, even a very large deviation of the assumed G value from the actual 
one would result in a negligible error. 

The data are well described with a Fourier polynomial of 2nd order. The 
resulting best-fit synodic period is Psyn = 103.02 ± 0.40 h. Pgyn agrees with 
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the Berthier et al. (2007) result of 102.94±0.11 h for the period of the mutual 
orbit, i.e. the rotational variability and the mutual orbit are synchronized. See 
Fig. 1 for a composite lightcurve in which the single-night observations are 
folded with the lightcurve period, along with the Fourier fit. The lightcurve 
amplitude is found to be 0.070 ± 0.005 mag peak-to-peak. 

To our knowledge, no complete lightcurve of Patroclus has been published 
previously, although it was known to have a long period and low amplitude 



(Angeli et al. , 1999). No signs of multiple periodicity (which might be ex- 
pected in a binary system) are discernible in our data. This implies that no 
mutual events occurred at the time of our observations (consistent with the 



Berthier et al. , 2007, model) and indicates a fully synchronized binary sys- 
tem. The small lightcurve amplitude implies that both components are nearly 
spherical in shape. While the data do not allow an unambiguous shape model 
to be determined, the observed lightcurve can be reproduced reasonably well 
assuming two slightly prolate spheroids with axial ratios 1.07:1:1, a realistic 



photometric function, and the Berthier et al. spin axis (see Fig. 1) 



3 Spitzer observations 



Patroclus was observed using the InfraRed Spectrograph (IRS; Houck et al 



2004) on board the Spitzer Space Telescope (Werner et al. , 2004). IRS was used 



in low-resolution spectroscopy mode using the modules SLl, LLl, and LL2 to 
obtain flux-calibrated spectra in the nominal wavelength range 7.4-38 pm at 
a relative spectral resolution A/AA between 64 and 128. The observed flux 
is practically purely thermal; the reflected component is negligible at these 



wavelengths (Emery et al. , 2006). The angular separation of Patroclus and 
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Table 2 

Observing geometry at the epoch of our observations. The hehocentric and Spitzer- 
centric coordinates are longitude and latitude in J2000 ecliptic coordinates. All 
values are constant during our observations to ±1 in the last quoted digit or better. 



The absolute visible magnitude equals H = 8.19 (quoted after Tedesco et al. , 2002 ); 



the slope parameter of the phase curve is assumed to be G = 0.15. 



Event: 


1 


2 


heliocentric distance r 


5.947 AU 


5.947 AU 


Spitzer-centric distance A 


5.95 AU 


5.98 AU 


solar phase angle a 


9.80° 


9.77° 


heliocentric coords 


170.8°, +18.03° 


170.9°, +18.00° 


Spitzer-centric coords 


160.5°, +18.2° 


160.7°, +18.1° 



Menoetius does not exceed 0.2" while the IRS pixel scale is 1.8" or coarser; 
the system was not spatially resolved. 

Time-resolved observations were obtained during two consecutive mutual events 
in June 2006, referred to as events 1 and 2 in the following. In event 1, Patro- 
clus shadowed Menoetius, and vice versa in event 2. The diameter ratio is only 
~ 1.1 and the mutual orbit is circular, hence the two events produced very sim- 
ilar observable effects. Both events lasted about 4 h and were pure shadowing 
events (the line of sight toward Spitzer was never obstructed). Thermal emis- 
sion from the shadowed surface areas therefore contributes to the observable 
flux. 



A total of 18 thermal-infrared spectra of Patroclus were obtained, nine per 
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Table 3 

Start times of our 18 Spitzer observations. See Fig. 2 for the respective system 
geometry. 





Day 


Time 




Day 


Time 




(June 2006) 


(UT) 




(June 2006) 


(UT) 


1.0 


24 


18:40 


2.0 


26 


10:42 


1.1 


24 


21:54 


2.1 


26 


23:22 


1.2 


24 


22:47 


2.2 


27 


00:24 


1.3 


24 


23:54 


2.3 


27 


01:31 


1.4 


25 


00:41 


2.4 


27 


02:19 


1.5 


25 


01:47 


2.5 


27 


03:29 


1.6 


25 


02:49 


2.6 


27 


04:24 


1.7 


25 


04:12 


2.7 


27 


05:55 


1.8 


25 


05:24 


2.8 


27 


06:52 



event, referred to as 1.0-1.8 and 2.0-2.8, respectively. See table 2 for the ob- 
serving geometry during the events. Observation start times are given in table 
3; each observation lasted ~ 6 min. As can be seen in Fig. 2, our observa- 
tions provide good temporal coverage of both events and their aftermath. The 



Marchis et al. (2006a) orbit model was used to predict the time of the events 



several months ahead. Observations 1.0 and 2.0 were performed well before 
the start of the events, in order to establish a flux baseline. 
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June 24 21:54(1.1) 

^» 


June24 22:47 (U) 


June 24 23:54 (1.3) 

• 


June 25 00:41 (1.4) 


4 






June 25 01:47 (1.5) 

• 
• 


June 25 02:49 (1.6) 

• 
• 


June 25 04:12(1.7) 

• 
• 


June 25 0524(1.8) 

• 

• 




June 26 23 22 (2.1) 

• • 


June 27 00:24 (22) 

• • 


June 27 01:31 (2.3) 


June 27 02:19 (2.4) 


June 27 03:29 (2.5) 

• 


June 27 04:24 (2.6) 

t 


June 27 05:55(2.7) 

• 
• 


June 27 06:52 (2.8) 

• 
• 



Fig. 2. Geometry of the system's mutual orbit during our Spitzer observations as 



seen from Spitzer, generated using the Berthier et al. (20071 orbit model. Each frame 



corresponds to one of our observations 1.1-1.8 (top two rows) and 2.1-2.8 (lower 
two rows), see table 3. Observations 1.0 and 2.0 are not plotted; at this scale, only 
one component would be displayed. Each frame is centered on Patroclus; its angular 
diameter is ~ 12.4 milli arcsec. During event 1, Menoetius is shadowed by Patroclus; 
and vice versa during event 2. The line of sight toward Spitzer was never obstructed. 



The event timing was shifted somewhat relative to the Berthier et al. prediction 
(within the uncertainties) to match the best-fit eclipse timing found in sect. 7.2. 
The system is eclipsed during observations 1.1-1.4 and 2.3-2.6, respectively. 

To prevent flux from 'spilling over' the edge of the slit, IRS spectroscopy 
targets must be accurately centered into the slit. The projected width of the 
LL slit is above 10", while that of the SLl slit is 3.7", comparable to the width 
of the point-spread function (PSF) in the respective wavelength range. For 
comparison: the pixel scale is 1.8 pm in the SL modules and 5.1 pm in the 
LL modules. "Blind" telescope pointings have a l(T-accuracy of 0.5" and are 
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therefore adequate to center sources of well known position (such as Patroclus) 
in the wide LL slits but risk placing a significant fraction of the PSF outside the 
SLl slit. Instead, small "spectral maps" were created for the SLl observations, 
with three small steps perpendicular to the slit, each offset by 2" (roughly half 
the slit width). This allowed us to estimate the target offset from the slit 
center and to correct for the effects thereof in the data analysis (see below). 

The usual "nod" strategy was used, i.e. the target was placed at about 33 and 
66 % of the slit length in each module in order to enable subtractive correction 
for diffuse background emission ("sky background"). 



The data obtained were reduced using the methods described in Emery et 



al. (2006, sect. 3). Briefly, the Spitzer Science Center (SSC) receives the data 
from the telescope and runs it through a reduction pipeline that converts 
raw data into 2-D spectral flux images (electrons/s) while flagging bad pixels, 
correcting known stray light artifacts, and performing fiat field corrections. 
The output is termed Basic Calibrated Data (BCD), which is described in 



detail in the IRS Data Handbook (http://ssc.spitzer.caltech.eciu/irs/ 



dh/). The data presented in this paper result from the BCD pipeline version 



S16.1.0. 

The BCD files contain background emission, mostly from the zodiacal cloud, 
which is removed by subtracting the two nod positions from one another for 
each observation. Extraction of the data to 1-D spectra essentially follows 
that of ground-based spectral data, where the object is identified in each row 
and the dispersion mapped across the chip. The signal (e~/s at this point) is 
the sum of the pixels within the chosen extraction width at each wavelength. 
Fractional pixels are included assuming the fiux is evenly distributed across 
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the individual pixel. The SSC provides a wavelength map of the chip for each 
spectral order. Since Spitzer is diffraction limited, the PSF, and therefore op- 
timal extraction width, varies with wavelength. We use the same extraction 
width as recommended by the SSC: ~ 2 x FWHM (14.4", 21.68", and 36.58" at 
12 pm, 16 pm, and 27 pm, respectively) and varying linearly with wavelength. 
Because of the varying PSF width and several other artifacts, the extracted 
spectral shape for each order is known to be slightly incorrect. The SSC cor- 
rects ("tunes") the final spectral shape using a 5th order polynomial fit to 
the ratio of stellar calibrator observations to models of their spectral shape 



(Decin et al. , 2004). These ratios are also used for final flux calibration (i.e., 
converting from e~/s to Jy). We perform basically the same procedure, with 
the exception that we use different calibrator stars, and we do not fit a poly- 
nomial, but rather apply the average ratio of the stars to the model directly. 
For the data presented here, we used del UMi, HR 2194, and HR 7891 for 
the SL modules and HD173511, del UMi, and HR 7341 for the LL modules. 
The different orders (SLl, LL2, and LLl) overlap slightly, and we scale the 
different orders relative to each other using these overlap regions. We then 
rescale the entire spectrum to the mean of the individual scale factors. The 
different low resolution orders generally agree within < 7 %. 



Observations at some wavelengths are compromised by defective detector pix- 
els and are disregarded in the following, e.g. all wavelengths above 33 pm. 
Systematic flux uncertainties due to, e.g., the residual effect of source mis- 
centering or uncertainties in the absolute flux calibration are estimated to 
be 3 % and added to the statistical uncertainty. The total flux uncertainty 
is dominated by the systematic uncertainty except for the low fluxes on the 
Wien slope at the shortest wavelengths. Final fluxes are provided in the On- 
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Fig. 3. Top: Ratio of two measured thermal spectra relative to pre-event measure- 
ment 1.0. Observation 1.2 (bottom) was taken during the eclipse, observation 1.7 
(top) in the aftermath. Note the difference in flux level and spectral slope, indica- 
tive of the eclipse-induced temperature drop. Bottom: Spectrum 1.0 in absolute flux 
units for comparison. 



line Supplementary Material. See Fig. 3 for a plot of representative flux ratios 
during event 1. 
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4 Emissivity features 



The observed spectra contain slight spectral features due to silicates within the 
wavelength ranges 10-12 pm and 18-22 pm. Emissivity spectra are calculated 
by dividing the measured flux spectrum by the modeled thermal continuum. 
Since we are not interested in smooth spectral slopes but in discrete features, 
the specific choice of model continuum is uncritical at this point. We see no 
reliable spectral variation among the 18 observations and therefore show the 
grand average. 

Our final emissivity spectrum result for the Patroclus system is shown in Fig. 4 
and compared to that of the Trojan asteroid (624) Hektor ( Emery et al.[ 2006[ ). 



Features near 10 and 20 pm are clearly visible; small scale (single channel) 
structure in the Patroclus spectra, e.g. near 25 pm, are residual calibration 
artifacts. 



Emery et al. conclude that the features near 10 and 20 pm in Hektor (and 
two other Trojan asteroids) are due to the presence of fine-grained (< few 
pm) silicates, likely in either a very under-dense ("fairy-castle") structure or 
embedded in a matrix material that is fairly transparent at these wavelengths 
(e.g., certain organic materials). Since the broad emissivity peaks in the Pa- 
troclus spectrum occur at the same wavelengths and with very similar shapes 
to those in the spectrum of Hektor, it is reasonable to infer similar mineralogy 
for the surface of Patroclus. The features, particularly the one near 10 pm, 
are more subdued in the spectrum of Patroclus than that of Hektor (or the 



other two Trojans presented in Emery et al. , 2006). This is most likely due 



to differences in grain size or packing state of the surface, two factors which 
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15 20 25 30 
Wavelength (|im) 



Fig. 4. Emissivity spectra of Patroclus (see text) and Hektor (Emery et al. , 2006 1 . 
Spectra are shifted vertically for clarity. 

strongly affect tlie sliape of spectra at tliese wavelengtlis, but could also imply 
different mixing ratios of silicates to opaque material. 



In order to avoid biases in thermal modeling, the wavelengths at which these 
features occur (10-12 pm and 18-22 pm) were disregarded in the thermal anal- 
ysis, leaving 178 data points per observation. 



5 Thermophysical modeling 

Analyzing thermal-infrared observations of asteroids requires a thermal model 
to calculate expected fluxes. To this end, the temperature distribution is calcu- 
lated based on the observing geometry and assumed physical properties, then 
the observable model flux is calculated by integrating the Planck function over 
the visible surface. 

Thermal fluxes depend on a number of physical properties, among them the 
effective diameter D, geometric albedo shape, thermal inertia F, spin pe- 
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riod P, spin-axis orientation, and amount of surface roughness. Typically, the 
absolute visible magnitude H is known, such that pv and D are linked through 
(see 



Pravec and Harris 2007 Appendix A, and references therein) 



D = 10-^^/5 1329 km ^ 



reducing the number of free parameters by one. 

The surface temperature distribution is governed by V (together with the spin 
rate): An object with zero thermal inertia (or vanishing spin rate) would dis- 
play a prominent temperature difference between the hot sub-solar point and 
the T = night side, whereas increasing F (or spin rate) shifts the tempera- 
ture maximum toward the afternoon side and, for observations at low phase 
angles, reduces the apparent color temperature, since solar energy is stored 
and re-emitted at night. 

Surface roughness, e.g. due to cratering, causes thermal-infrared beaming: 
relative to a smooth surface, a rough surface displays more surface elements 
facing the Sun, even close to the limb, leading to a preferentially sunward 
emission and, for observations at low phase angle, a higher color temperature. 
Also, surface elements inside an indentation can radiatively exchange energy, 
leading to mutual heating. 

Simple thermal models in which highly idealized assumptions are made on 
shape, thermal inertia, and surface roughness have enjoyed widespread us- 
age over the past decades and have been shown to allow robust estimates of 



diameter and albedo in many cases (see Harris , 2006 , for a recent overview) 



Deriving the thermal inertia, on the other hand, requires more realistic ther- 
mophysical modeling, and in particular some knowledge of the object's spin 
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state, which is currently available for a small (but growing) number of asteroids 
only. Also, given the notorious difficulty of ground-based mid-IR observations, 
complex models with a large number of free parameters are not frequently 
warranted by the data quality. Nevertheless, detailed thermophysical models 
(TPMs) in which shape, spin state, thermal inertia, and surface roughness 



are explicitly taken into account, have been developed (e.g. Spencer, 1990 



Lagerros 


1996 


1998 


Mueller 


2007 



used successfully to analyze mid-IR observations of asteroids (e.g. Miiller and 



Lagerros 


1998 


Harris et al. 


2005 


2007 


Mueller 


2007) 



Since the thermal skin depth, a measure of the penetration depth of the heat 
wave, is typically in the cm-range, lateral heat conduction can be neglected. We 
also assume, as is common practice, that all relevant parameters are constant 
with depth and temperature^ It is then practical to use dimensionless units 
(see 



Spencer et al. , 1989) for time t', depth X, and temperature T'. The only 



free parameter in the heat conduction problem is then the thermal parameter 



6, as defined by [Spencer et al.[ which is proportional to T. The equations 
governing thermal conduction into the subsoil are 



(X, t') = (X, t') 

dt' ^ ' ' (9X2 V , ; 

T'{Q,t'f = ^^s{t') + Q^T'M 



(2) 
(3) 
(4) 



^ In the case of lunar regolith, thermal parameters are actually empirically known 



to depend on both temperature and depth (see, e.g., Jones et al. 1975 , and references 
therein). It is, however, unclear how to apply lunar results to asteroid regolith, which 
may be substantially different. Thermal- inertia results derived from TPMs should 
be thought of as effective averages over the relevant length scales. 
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where ns is the cosine of the local solar zenith distance (clipped to be > 0). 
Eqn. 2 describes heat diffusion within the subsoil, eqn. 3 the surface boundary- 
condition (thermally emitted energy equals absorbed solar energy plus the net 
heat flux from the subsoil), eqn. 4 ensures the regularity of the solution. See, 



e.g., Mueller (2007, Sect. 3.2.2) for a more detailed discussion. 



Surface roughness can be modeled (following Spencer et al. , 1989 Spencer 



1990) by multiplying surface temperatures as calculated from Eqns. 2-4 with 



a factor rj' (//' is called (3 in the quoted papers), rj' < 1 leads to elevated 
apparent color temperature. Increasing roughness leads to decreasing values 



of r/', a smooth surface would display r/' = 1. Spencer et al. (1989) derive 
Tj' ~ 0.72 for the lunar surface. 



5.1 Binary TPM (BTPM) 



We describe a new Binary TPM (BTPM). It is based on the Spencer (1990) 
TPM, but includes the effect of mutual events (eclipses and occultations) on 
thermal fluxes. Our BTPM is limited to fully synchronized binary systems 
with a circular mutual orbit, such as Patroclus-Menoetius. Such objects are 
at rest in a co-rotating coordinate system. We can therefore use a single-body 
TPM to model the binary system, using a shape model which consists of two 
disjoint parts. 

In contrast to most TPMs, which assume a globally convex shape, surface 
elements can now shadow one another and obstruct the line of sight toward 
the observer. In principle, some facets can also radiatively exchange energy, 
similar to mutual heating inside craters. That latter effect is not included in the 
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model. While it is important in the thermal modeling of globally non-convex 
bodies, it is negligible for eclipse observations at low solar phase angle. 



When generating the binary shape model, both components are modeled as 
a mesh of triangular facets. Arbitrary component shapes are allowed by the 
model code; spheres and ellipsoids are used for this study. Each facet is checked 
for potential shadowers, i.e. other facets which appear above the local hori- 
zon. The list of potential shadowers along with details on the mutual viewing 



geometry is stored as part of the shape model (see Mueller, 2007, Sect. A. 1.1 
for details), ensuring that these calculations do not need to be repeated for 
each BTPM run. 

Since mutual heating is neglected, the observable emission can be calculated 
facet-by-facet. Mutual heating would introduce an interdependence of surface 
temperatures and add to the complexity. Temperatures are calculated by nu- 
merically solving Eqns. 2-4. Eclipses are taken into account by setting fis, the 
insolation term in eqn. 3, to zero when the facet in question is eclipsed by 
another facet. Facets are considered eclipsed when their midpoint is shadowed 
(no partial shadowing of individual facets). Care must be taken to choose a 
sufficient time resolution, so as to provide appropriate resolution of eclipse 
events, which are short compared to the rotation period. 



6 NEATM analysis 



As a first step, the obtained spectra (minus the location of emissivity fea- 



tures) were analyzed using the NEATM (Harris, 1998). For the purposes of 
this study, it is useful to think of the NEATM as a simplified TPM in which 
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the shape is assumed to be spherical and thermal inertia is not directly mod- 
eled (therefore the spin state is irrelevant). It contains a model parameter 77, 
which modifies model surface temperatures in the same way as our BTPM's v( 
(temperature oc 77"^/^), but is used for an effective description of both surface 
roughness and thermal inertia. For observations at low phase angles, increas- 
ing thermal inertia lowers surface temperatures and therefore increases r] (all 
other parameters kept constant). 

The NEATM has been used widely to analyze thermal-IR observations of 
asteroids and is known to provide diameter results which are in generally 
good agreement with diameter determinations obtained using other methods. 
Additionally, conclusions on thermal properties such as thermal inertia can be 



drawn from the best-fit value. See, e.g., Harris (2006) for a recent overview. 



We fitted the NEATM to each of our 18 observations, see table 4 for the 
results. Since the NEATM does not include the effects of eclipses, it is not 
applicable during and just after eclipses. By averaging over observations 1.0, 
1.7, 1.8, and 2.0-2.2 (see Fig. 2 and table 3), we obtain D ~ 146 km and 
?7 ~ 0.85. The 77 value indicates a relatively low thermal inertia and significant 
surface roughness, but does not provide a quantitative constraint on either. 

While the NEATM takes no direct account of shadowing, the reduced total 
fiux and the reduced temperature during and just after the eclipse events 
should be refiected in the NEATM results. Indeed, as can be seen in Fig. 5, 
there is a clearly recognizable dip in NEATM diameter during event 1 (results 
for event 2 are qualitatively identical), which refiects the event-induced fiux 
drop. Simultaneously, the best-fit 77 rises, corresponding to a lower apparent 
color temperature due to the echpse-induced cooling of the shadowed parts. 
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Table 4 

NEATM fit to our Spitzer data: best-fit r), diameter D, and geometric albedo py. 
Uncertainties were derived through a straight-forward Monte-Carlo analysis (mean 
and average of fits to 300 random spectra per observations, normally distributed 
about the measured data) and account for the statistical uncertainty, only. 



Obs. 


V 


D{km) 




Obs. 


rj 


D{km) 




1.0 


0.845 ± 0.004 


146.8 ±0.5 


4.34 ±0.03 


2.0 


0.839 ±0.004 


148.5 ±0.5 


4.24 ±0.03 


1.1 


0.852 ± 0.004 


143.5 ±0.6 


4.54 ±0.04 


2.1 


0.856 ± 0.004 


144.9 ± 0.6 


4.45 ± 0.04 


1.2 


0.864 ± 0.004 


141.3 ±0.5 


4.68 ±0.04 


2.2 


0.855 ±0.005 


143.4 ±0.6 


4.55 ± 0.04 


1.3 


0.860 ± 0.006 


140.9 ±0.7 


4.71 ±0.05 


2.3 


0.878 ±0.005 


141.5 ±0.5 


4.67 ± 0.04 


1.4 


0.856 ± 0.004 


143.2 ±0.5 


4.56 ± 0.04 


2.4 


0.895 ±0.006 


138.9 ±0.6 


4.85 ± 0.04 


1.5 


0.848 ± 0.004 


146.7 ±0.6 


4.34 ± 0.04 


2.5 


0.863 ± 0.004 


133.6 ±0.4 


5.24 ± 0.04 


1.6 


0.827 ± 0.004 


144.1 ±0.6 


4.50 ± 0.04 


2.6 


0.891 ± 0.004 


142.5 ±0.5 


4.60 ± 0.04 


1.7 


0.835 ± 0.004 


146.1 ±0.5 


4.38 ±0.03 


2.7 


0.875 ± 0.004 


144.1 ±0.6 


4.50 ± 0.04 


1.8 


0.851 ±0.004 


147.2 ±0.5 


4.32 ±0.03 


2.8 


0.861 ± 0.004 


143.3 ±0.6 


4.55 ± 0.04 



While no quantitative conclusion on thermal inertia can be drawn from the 
NEATM analysis, our results indicate that we have indeed observed the ther- 
mal response to eclipse events. 
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Fig. 5. NEATM fits to observations 1.0-1.8 (see table 4). Above: Best-fit effective 
diameter against observing time. Below: Best-fit r) against observing time. The 
vertical lines denote beginning and end of the eclipse event. The error bars are 
purely statistical. 

7 BTPM analysis 

7.1 Fitting technique 

In this section, we describe fits to our Spitzer data (denoted as di for the data 
and (Tj for their uncertainties) using the new BTPM described in sect. 5.1. 

Variable parameters are T, rj', area-equivalent system diameter Dj^, and the 
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time of the events (eclipse time predictions from the orbit model are currently 
uncertain by ± a few hours). To this end, synthetic lightcurves were gener- 
ated using the observing geometry given in table 2 (separately for each event) 
and for a large range of F and v( values at those wavelengths Aj at which 
Spitzer data have been obtained. The positions of emissivity features (10-12 
and 18-22 pm, see sect. 4) were disregarded. Model fluxes mj(Ai) were calcu- 
lated for 1000 time points per revolution, corresponding to a time resolution 
of ~ 6 min or roughly the length of one Spitzer observation. This provides 
more than 40 time points across the duration of each ~ 4 h event. The rough 
rotational phase of the events was determined from visual inspection of the 
synthetic lightcurves, the exact timing was determined by fitting a time offset 
At. All model calculations were performed assuming H = 8.19, = 0.0433, 
and = Do = 147 km. Fluxes are assumed to be proportional to Dj^: 
mi{Xi,DX) = X mj(Aj,Do) with x = (D_4/Do)^j3 For each considered set of 
F, 1]' , and At the best-fit x is found by minimizing (through linear regression) 

The corresponding best-fit diameter equals = ^/x 147 km. 

After initial runs over wider and coarser grids, an equidistant grid was used 
with a step width of 0.5 J s'^^^K'^m"^ in F, 0.0025 in rj', and 1 min in At. Data 
were analyzed separately for each event. For each combination of F, rj', and 



^ Implicitly, we neglect the albedo dependence of thermal fluxes: varying Z?_4 with 
constant H varies py. Temperatures depend on the absorptivity, which is a strong 
function of py for high albedo, but a weak function for low py. Since Patroclus' 



albedo is very low, we reckon this approximation is uncritical (see also Harris and 



Harris, 1997 1 . 
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At, the best-fit x and corresponding were determined and stored. The grid 
was then searched for that combination of F, 77', and At which led to the global 
minimum in x^- 



In order to study the accuracy of the results, a Monte-Carlo technique was em- 
ployed. For each observation, 5000 random ("noisy") spectra were generated; 
for each wavelength the distribution of the 5000 random points is Gaussian, 
with a mean value equal to the measured flux and a standard deviation equal 
to the flux uncertainty. Best-fit BTPM parameters were determined for each 
set of nine random spectra; mean value and standard deviation of the resulting 
values were adopted as nominal value and statistical uncertainty, respectively. 
This error analysis does not yet account for the uncertainty in the axis of the 



mutual orbit (A = 241.37 ± 0.33°, p = -57.35 ± 0.36°; [Berthier et al] |2007| ). 
Eclipse depth depends sensitively on the sub-solar latitude and therefore on 
the orbital axis, thus the axis uncertainty (while small) could be expected to 
have a significant impact on the accuracy of our thermal-inertia result. We re- 
peated the Monte-Carlo BTPM analysis assuming the axis within the quoted 
range of uncertainty that leads to the largest change in sub-solar latitude: 
A = 241.70°, j3 = -57.71°. The offset axis leads to a somewhat increased sub- 
solar latitude, implying shallower eclipse events, i.e. a smaller area fraction of 
the eclipsed component is shadowed. 



Additional simulations were performed in which the components were assumed 
to be prolate spheroids with aligned long axes (nominal orbit solution). Axis 
ratios (identical for the two components) of up to 1.09 were tried; larger ratios 
would be inconsistent with the low visible lightcurve amphtude (see sect. 2). 
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Table 5 

Best-fit BTPM parameters for events 1 and 2. The two lines of results per event 
refer to the nominal orbit solution and that offset by la (see text), respectively. 
Syzygy times are as seen on Spitzer. There are 9 x 178 = 1602 data points per event 
and four fit parameters, hence a best- fit iii the low thousands is to be expected. 



Syzygy 



{JD - 2453910) 



X 



(J s-i/2K-im-2) 



Da 



(km) 



Event 1 



Event 2 



1.466 ±0.003 2543 ±84 



1.470 ±0.002 2586 ±; 



3.624 ± 0.001 3566 ± 102 



3.624 ±0.001 3847 ±109 



20.7 ±3.8 0.762 ±0.014 145.7 ±0.3 



7.6 ±1.7 0.814 ±0.008 146.0 ± 0.3 



6.4 ±0.9 0.838 ±0.005 143.4 ± 0.3 



5.1 ±0.8 0.845 ±0.004 143.1 ±0.3 



7.2 Results 

BTPM results are reported in table 5 for both the nominal and the offset axis. 
Results from the ellipsoidal-shape model runs are indistinguishable and are 
not reported. 

The best-fit diameter values resulting from the four BTPM analyses and the 
NEATM analysis (~ 146 km; see sect. 6) are mutually consistent at the 1 % 
level and average to 144.6 km. 



The time of the syzygy (maximum eclipse depth, when the two components 
and the Sun are aligned) was fitted to the data and is consistent with the 



Berthier et al. (2007) model prediction within the uncertainty of ± a few 



hours. The times given in table 5 are for a Spitzer-centric system; subtract 
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the one-way light-time for the corresponding times at Patroclus. 



Model runs with the offset axis lead to lower thermal-inertia results in both 
cases. This is expected since the offset reduces the eclipse depth: The leading- 
order effect of increasing thermal inertia or decreasing eclipse depth is to 
reduce the echpse-induced flux drop; hence underestimated eclipse depth leads 
to underestimation of thermal inertia. Final error estimates are obtained by 
adding in quadrature the statistical uncertainty in the nominal solution to the 
difference between the nominal and offset-axis solutions; our final results are 
r = 21 ± 14 J s^/^K-^m^ (event 1) and T = 6.4 ± 1.6 J s-^^^K-^m^ (event 2). 

While our two thermal-inertia results differ by a factor > 3, they are mu- 
tually consistent given the uncertainties. Note, however, that neither result 
is representative of the system as a whole, but is dominated by the respec- 
tive shadowed region. A difference between the two thermal-inertia results 
would therefore imply inhomogeneous regolith properties and might hint at 
component-to-component differences. While tantalizing, our data do not al- 
low this conclusion to be drawn at a statistically significant level. We estimate 
the system-average thermal inertia to be between 5 and 35 J s"^/^K"^m"^ or 
20 ± 15 J s'^/^K'^m"^, where the scatter may be partially due to surface inho- 
mogeneity. 

Our thermal-inertia results depend sensitively on the assumed spin axis, as 
evidenced by our experiment with the offset axis and, more drastically, by our 
experience with preliminary orbit models. E.g., in a preliminary data analysis 



based on a preliminary orbit model (Mueller, 2007, Sect. 6.8), a best-fit T of 



~ 90 J s"^/^K"^m"^ was obtained. The uncertainty introduced by the spin-axis 
uncertainty was not studied in the previous analysis. We point out, however. 
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that the new orbit model fits the Spitzer data much better: is now reduced 
by a factor > 3 and is now in the expected range given the number of data 
points. Also, the nominal spin axis fits the Spitzer data better than that offset 
by la, inspiring further trust in the presented results. 

The best-fit rj' ranges from 0.76 to 0.85, between the lunar value of rj' ~ 0.72 



(Spencer et al. , 1989) and the NEATM best-fit r] ~ 0.85, which describes the 
(cooling) effect of thermal inertia in addition to roughness. The best-fit rj' is 
inversely correlated with best-fit F. The impact of thermal inertia on surface 
temperatures is twofold: It changes the response to the eclipse (i.e. depth, 
shape, and wavelength dependence of the eclipse- induced flux drop), but it 
also determines the diurnal temperature distribution and therefore the overall 
spectral shape (color temperature). As will be discussed in the next section, 
the effective thermal inertia for these two phenomena need not be the same 
because of the different time and depth scales involved. Our model, however, 
assumes a homogeneous thermal inertia. Due to the design of our observations, 
the data analysis is more sensitive to the eclipse effect than to the diurnal 
effect; we take the variation in best-fit r]' to be an effective description of the 
residual diurnal effect rather than an indication of roughness inhomogeneity. 



8 Discussion 



8.1 Thermal inertia 

Our thermal-inertia result is below the value for lunar regolith (50 J s"^/^K"^m"^), 
suggesting a very fine regolith and dearth of rocky outcrops. For a meaning- 
ful comparison, however, the lower temperatures at a heliocentric distance of 
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r ~ 6 AU must be taken into account: For predominantly radiative heat trans- 
fer, as expected in a fine regolith on an airless body, the thermal conductivity 
K scales with T^, hence 

r oc oc T^/^ oc r"^/^ (6) 

The thermal inertia of Patroclus (20 ± 15 J s'^^^K'^m"^ at 5.95 AU) scaled 
to r = 1 AU is around 76 J s'^^^K'^m"^, comparable to the lunar value and 
suggesting a similarly fine regolith. 

Our result is similar to published upper limits on the thermal inertia of two 



other Trojans (Fernandez et al. , 2003, see also table 6) which were, however. 



based on a significantly less extensive database and a highly indirect method. 
We are not aware of any other thermal-inertia determination of Trojans. 

Table 6 contains a list of published thermal-inertia values for Jovian satellites, 
which are at the same heliocentric distance as Patroclus. Two values are given 
for lo, Ganymede, and Callisto, respectively, based on ground-based eclipse 
observations (first line) or on thermal imaging of the diurnal temperature 
distribution from near-by spacecraft (second line). The latter are effective disk- 
average values; significant thermal-inertia inhomogeneity was found across the 
surfaces of the satellites. Patroclus' thermal inertia is roughly in line with the 
satellites' eclipse thermal inertia, which is lower than their diurnal thermal 
inertia by a factor of several. Different explanations for this difference in the 



satellites have been published; Greeley et al. (2004) review two mechanisms 



which require significant amounts of surface ice, and a third mechanism (first 



proposed by Morrison and Cruikshank, 1973) based on vertical thermal-inertia 



inhomogeneity: Because eclipse events are much shorter than the spin period, 
the corresponding thermal wave penetrates much less deep. Keeping in mind 
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Table 6 

Thermal-inertia measurements of Trojans (top 3 lines) and Jovian satellites (below). 



Thermal inertia Method Reference 



(J s-i/2K-im-2) 



(617) Patroclus 



(2363) Cebriones 



(3063) Makhaon 



Europa 



lo 



Ganymede 



Callisto 



20 ±15 



< 14 



< 30 



45-70 



13 ±4 



70 



14 ±2 



70 



11 ± 1 



50 



eclipse this work 



diurnal 



diurnal 



diurnal 



eclipse 



diurnal 



eclipse 



diurnal 



Fernandez et al. 



Fernandez et al. 



Spencer et al. 



(20031 



(20031 



(19991, 



Greeley et al. 



Morrison and Cruikshank 



diurnal Rathbun et al. 



(20041, 



(19731 



(20041 



McEwen et al. 



(20041 



eclipse Morrison and Cruikshank 



Spencer 



(1987) 



(19731 



Pappalardo et al. 



Morrison and Cruikshank 



Spencer 



(1987) 



(19731 



Moore et al. 



(20041 



(20041 



that the TPM 'thermal inertia' is really an average over the relevant depth 
scale (see footnote 1 on p. 24), the effective diurnal thermal inertia is averaged 
over deeper regions than its eclipse counterpart. The thermal-inertia difference 
then implies that near-surface material has a lower thermal inertia than the 
material underneath (at cm depths), consistent with expectations that the 
most under-dense and 'fluffiest' material is located on top of more compact 
material. Such vertical grain-size sorting has been observed directly in the 
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lunar regolith during the Apollo era. 



For Patroclus, we do not have a measurement of the diurnal thermal inertia, 
but the eclipse thermal inertia is similarly low. No evidence of surface ice 



has been found in the spectra of Patroclus (Emery and Brown, 2003, 2004), 
hence the two ice-based mechanisms to increase the diurnal thermal inertia 
do not apply. A measurement of the diurnal thermal inertia, combined with 
our result, would therefore allow depth-resolved information to be obtained to 
some extent. 



While comparisons with the Galilean satellites are instructive, we caution that 
Patroclus is smaller than the satellites by more than an order of magnitude. 
For an asteroid sample ranging from sub-km near-Earth asteroids to (1) Ceres, 



Delbo' et al. (2007) found a trend whereby larger objects have lower thermal 



inertia. 



Mueller et al. (2006) and Delbo' and Tanga (2009) report constraints 



on the thermal inertia of D ~ 100 km main-belt asteroids (Lutetia, Lydia, 
Ekard), similar to Patroclus in size. While the uncertainties are large, F ~ 
100 J s"^/^K"^m"^ appears to be a typical value in this size range. Using eqn. 6, 
this corresponds to F ~ 60 J s"^/^K"^m"^ at r ~ 6 AU, significantly above our 
result but consistent with the diurnal thermal inertia of the Galilean satellites. 
This is consistent with either a systematically finer regolith on Patroclus (it 
is unclear what could cause this) or the above-mentioned vertical grain-size 
sorting. 
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8.2 Diameter and mass density 



The diameter results from two independent data analyses (NEATM and BTPM) 
agree at the 1 % level and average around ~ 145 km for the area-equivalent 
diameter of the system as a whole: D^^ = Df + D|. The statistical diameter 
uncertainty is negligible relative to the systematic uncertainty, although the 
latter is hard to estimate. Since it is more realistic, the BTPM would be ex- 
pected to have lower systematics. While more work is required to accurately 
gauge the systematic uncertainty in TPM-derived diameters, it was found not 



to exceed 10 % in the case of near-Earth asteroids (see Mueller, 2007). Due to 
its more regular shape and the low phase angle of our observations, Patroclus 
is less challenging to model than near-Earth asteroids. We conclude that 10 % 
is a conservative upper limit on the diameter uncertainty, corresponding to 
20 % in albedo (see eqn. 1). Our final results are: = 145 ±15 km and 
= 0.045 ± 0.009. 



Our result is intermediate between previous estimates: Tedesco et al. (2002) 



obtained 140.9 ± 4.7 km (assuming rj = 0.756), and Fernandez et al. (2003) 
obtained 166.0 ± 4.8 km (assuming rj = 0.94). The quoted error bars are 
purely statistical. We point out that their results were subject to significant 
systematic uncertainties due to the assumed (not fitted to the data) rj values. 
Systematic uncertainties are acknowledged by [Fernandez et al." to limit the 



accuracy of their results, but no quantitative discussion is given. Since our 
database has a significantly larger wavelength coverage and since we use a more 
realistic thermal model, the true uncertainty in the cited diameters should be 
large compared to that in our results. With this in mind, all diameter results 
are in reasonable agreement with one another. 
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Assuming a diameter ratio of D1/D2 = 1.082, the components' diameters 
are Di = 106 ± 11 and D2 = 98 ± 10 km. The corresponding bulk density is 
1.08 ±0.33 g cm"^, in excellent agreement with the Berthier et al. (2007) result 



of 0.90to;33 g cm"^ and indicative of a bulk composition dominated by water 
ice. Significant contributions by heavier materials (such as rock) are possible 
if compensated by interior voids. Near-infrared spectra show no indication 
of water ice on the surface, even in the 3 pm region, where even very small 



amounts of H2O would be detected (Emery and Brown, 2003, 2004). If the 



bulk composition does include significant H2O, the surface must be coated 
by a devolatilized mantle. This would be expected for a extinct-comet-like 
object, i.e. an icy object that spent significant amounts of time in the inner 
Solar System. 

It is not unusual for published asteroid diameters to be quoted with purely 
statistical uncertainties, neglecting the typically dominant systematics. Most 
known asteroid diameters are derived from thermal-infrared observations. While 
more research is needed to quantify their systematic diameter uncertainty, 
10 % seems to be an appropriate value, larger when the color temperature is 
assumed rather than fitted to the data, as is typically done in the widely used 
'Standard Thermal Model' (STM). This includes the entire SIMPS catalog 



(Tedesco et al. , 2002), i.e. the majority of currently known asteroid diameters. 



A realistic evaluation of the size uncertainty — including systematics — is crucial 
when published diameters are used to determine the mass density. Due to the 
dependence on D~^, density uncertainties are frequently dominated by that 
in diameter. Ignoring the systematic diameter uncertainty may lead to overly 
optimistic density estimates. A diameter uncertainty of 10 % translates into 
30 % in density; 15 % in diameter means some 45 % in density. 
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8. 3 Mineralogy 



It is instructive to compare Patroclus to its mythological antagonist, (624) 
Hektor, the only other currently known Trojan binary and the only other 
Trojan with known mass density. 

Patroclus and Hektor have nearly identical albedo and mid-IR emissivity spec- 
tra, yet their densities differ by a factor > 2: p ~ 2.2-2.5 g/cm^ for Hektor 



e.g. Weidenschilling , 1980 Marchis et al. , 2006b Lacerda and Jewitt 2007) 



to Patroclus' p ~ 1 g/cm^. Reflectance spectra of both are featureless from 
~ 0.3 to 4.0 pm, but Hektor displays a much steeper spectral slope throughout 
that wavelength range. It is uncertain whether the difference in spectral slope 
indicates distinct surface compositions or is a manifestation of differences in, 
e.g., grain size, abundance of opaque materials, or degree of space weathering. 



Trojans in general form two distinct groups based on visible and near-IR colors 



2007; 


Emery et al. 


2007; 


Melita et al. 


2008 


Roig 



et al. , 2008), one group similar to Hektor and one similar to Patroclus. It 
is tempting to take the density difference between Patroclus and Hektor as 
indicating a compositional difference between the two groups. In the light of 
this, their very similar mid-IR spectra are intriguing. 



The mutual orbits of the Patroclus and Hektor systems are conspicuously 
dissimilar: While Patroclus is a double-asteroid system with two components 
of roughly equal size and nearly spherical shape, Hektor consists of a large 
irregular primary (which is thought to be a contact binary itself) and a much 



smaller moonlet ( Marchis et al. 2006b ) . 
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9 Conclusions 



Patroclus Patroclus has a low thermal inertia of 20 ± 15 J s"^/^K"^m"^. This 
first thermal- inertia measurement for a Trojan indicates a surface covered in 
fine, mature regolith similar to large main-belt asteroids and the Galilean 
satellites. Independent evidence for a fine regolith is provided by the mid- 
IR emissivity features reported herein, which imply the presence of very fine 
silicate grains on the surface. 

Patroclus' optical lightcurve implies a nearly spherical shape for both com- 
ponents and indicates that the system is fully synchronized, probably due to 
tidal component- component interactions. 

We confirm that Patroclus' albedo is low (pv = 0.045 ± 0.009) and that the 
system's bulk mass density is 1.08 ± 0.33 g cm"^. Patroclus' known physical 
properties (albedo, density, fine silicate regolith, and published optical spectra 
which are similar to cometary surfaces but devoid of any traces of water ice) 
match those expected in extinct comets and suggest that Patroclus formed 
beyond the 'snow line' but spent significant amounts of time in the inner 
Solar System, leaving its surface devolatilized. Such a dynamical history is 



expected for Trojans in the framework of the 'Nice Model' (Morbidelli et al 



2005, 2009). While Patroclus is arguably the best-studied Trojan asteroid, it 
is unclear how representative it is for the Trojan population as a whole; the 
composition of the only other Trojan with known mass density, Hektor, is 
dominated by materials heavier than water ice. 
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Binciry eclipses in the thermal infrared We demonstrate the viability 
of a new method to measure the thermal inertia of asteroids, through thermal 
observations during and after shadowing events in binary systems. This new 
method is considerably more direct than the usual method of sampling the 
diurnal temperature distribution through observations at a large phase-angle 
range. 

In either case, the measured thermal inertia is an effective average over the 
relevant depth scale, which is smaller for eclipses than for rotation due to the 
shorter time scales involved. Vertical grain-size sorting may cause the eclipse 
thermal inertia to be lower than the diurnal thermal inertia. While this makes 
comparisons harder, it may shed light on the elusive vertical regolith structure. 

For objects in the outer Solar System, sufficient phase-angle coverage for the 
diurnal measurement method is hard or impossible to obtain. The maximum 
phase angle ccmax decreases with increasing heliocentric distance r (in AU): 
Q^max = arctan(l/r). For Neptune, e.g., ctmax ^ 2°. Thermal eclipse observa- 
tions may be the only method to measure the thermal inertia of Centaurs and 
Kuiper Belt objects from (close to) Earth - once the mutual orbits of such 
binary systems are reliably known. 
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